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1. Introduction 

Time and phase are not well defined concepts at the quantum level. The present work 
belongs to a longstanding effort to model phase noise and phase-locking effects that 
are found in highly stable oscillators. It was unambiguously demonstrated that the 
observed variability, (i.e., the f/f frequency noise of such oscillators) is related to the 
finite dynamics of states during the measurement process and to the precise filtering 
rules that involve continued fraction expansions, prime number decomposition, and 
hyperbolic geometry 

We add here a quantum counterpart to these effects by studying the notions of 
quantum phase-locking and quantum phase entanglement. The problem of defining 
quantum phase operators was initiated by Dirac in 1927 [5 j. For excellent reviews, see 

We use the Pegg and Barnett quantum phase formalism 7j where the calculations 
are performed in an Hilbert space H q of finite dimension q. The phase states are defined 
as superpositions of number states from the so-called quantum Fourier transform (or 
QFT) 

|fp> = g- 1/a Eexp(^)|n>. (1) 

n=0 " 

in which i 2 = —1. The states \9 P ) form an orthonormal set and in addition the projector 
over the subspace of phase states is J2l=o \@p) {Qp\ — lg where l q is the identity operator in 
H q . The inverse quantum Fourier transform follows as |n) = g -1 ^ 2 J2p=o^ x p{— 2lnpn )\9 p ). 
As the set of number states \n), the set of phase states \9 P ) is a complete set spanning 
H q . In addition the QFT operator is a q by q unitary matrix with matrix elements 

From now we emphasize phase states \9') satisfying phase- locking properties. The 
quantum phase-locking operator is defined in ([19)1 . We first impose the coprimality 
condition 

(p,?) = l, (2) 

where (p, q) is the greatest common divisor of p and q. Differently from the phase states 
the \9' p ) form an orthonormal base of a Hilbert space whose dimension is lower, 
and equals the number of irreducible fractions p/q, which is given by the Euler totient 
function 0(g). These states were studied in our recent publication 

Guided by the analogy with the classical situation [2|, we call these irreducible 
states the phase-locked quantum states. They generate a cyclotomic lattice L [E] 
with generator matrix M of matrix elements k'S$, (p,q)=l and of size <p(q). The 
corresponding Gram matrix H = M^M shows matrix elements h^] = c q (n — I) which 

|| Some errors or misunderstandings are present in that earlier report. The summation in (3), (5), (7) 
and (9) should be (this is implicit) from to <fi(q). The expectation value (0 l ° ck ) in (8) should be 
squared. There are also slight changes in the plots. 
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are Ramanujan sums 



c q (n) = exp(2i7r-n) 



p 



<f>{<h) 



with qi = q/ (q, n). 



(3) 



where the index p means summation from to q — 1, and (p, g) = 1. Ramanujan sums 
are thus defined as the sums over the primitive characters exp(227r— ), (p,q)=l, of the 
group Z q = Z/qZ. In the equation above fi(q) is the Mobius function, which is if the 
prime number decomposition of q contains a square, 1 if q = 1, and (— l) k if q is the 
product of k distinct primes [§]. Ramanujan sums are relative integers which are quasi- 
periodic versus n with quasi-period 0(g) and aperiodic versus q with a type of variability 
imposed by the Mobius function. Ramanujan sums were introduced by Ramanujan in 
the context of Goldbach conjecture 0. 

They are also useful in the context of signal processing as an arithmetical alternative 
to the discrete Fourier transform ^U]. In the discrete Fourier transform the signal 
processing is performed by using all roots of unity of the form exp(2iirp/q) with p 
from 1 to q and taking their nth powers e p (n) as basis function. We generalized the 
classical Fourier analysis by using Ramanujan sums c q (n) as in (JSJ) instead of e p (n). 
This type of signal processing is more appropriate for arithmetical functions than is the 
ordinary discrete Fourier transform, while still preserving the metric and orthogonal 
properties of the latter. Notable results relating arithmetical functions to each other 
can be obtained using Ramanujan sums expansion while the discrete Fourier transform 
would show instead the low frequency tails in the power spectrum. 

In this paper we are also interested in pairs of phase-locked states which satisfy the 
two conditions 



Whenever it exists p is uniquely defined from minus the inverse of p modulo q. 
Geometrically the two fractions p/q and p/q are the ones selected from the partition 
of the half plane by Ford circles. Ford circles are defined as the set of the images of 
the horizontal line z=x+i, x real, under all modular transformations in the group of 
2x2 matrices SL(2, Z) 4]. Ford circles are tangent to the real axis at a Farey fraction 
p/q, and are tangent to each other. They have been introduced by Rademacher as an 
optimum integration path to compute the number of partitions by means of the so-called 
Ramanujan's circle method. In that method two circles of indices p/q and p/q, of the 
same radius are dual to each other on the integration path [see also Part 2 and Fig. 



(j), q) = 1 and pp 



l(mod q). 



(4) 



3]. 



2. Hints to the hyperbolic geometry of phase noise 



2.1. The phase-locked loop, low pass filtering and 1/ / noise 



A newly discovered clue for 1/f noise was found from the concept of a phase locked 
loop (or PLL) [2]. In essence two interacting oscillators, whatever their origin, attempt 
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to cooperate by locking their frequency and their phase. They can do it by exchanging 
continuously tiny amounts of energy, so that both the coupling coefficient and the beat 
frequency should fluctuate in time around their average value. Correlations between 
amplitude and frequency noises were observed jTTj. 

One can get a good level of understanding of phase locking by considering the case 
of quartz crystal oscillators used in high frequency synthesizers, ultrastable clocks and 
communication engineering (e.g., mobile phones). The PLL used in a FM radio receiver 
is a genuine generator of 1/f noise. Close to phase locking the level of 1/f noise scales 
approximately as a 2 , where a = oK/iob is the ratio between the open loop gain K 
and the beat frequency uJb times a coefficient a whose origin has to be explained. The 
relation above is explained from a simple non linear model of the PLL known as Adler's 
equation 

9(t) + KH{P) sin 0(t) = lob, (5) 

where at this stage H(P) = 1, ujb = w{t) — ojq is the angular frequency shift between 
the two quartz oscillators at the input of the non linear device (a Schottky diode 
double balanced mixer), and 9{t) is the phase shift of the two oscillators versus time t. 
Solving (jSJ) and differentiating one gets the observed noise level a versus the bare one 
a = 5ujb/lJb- Thus the model doesn't explain the existence of 1/f noise but correctly 
predicts its dependence on the physical parameters of the loop [2j. 

Besides one can get detailed knowledge of harmonic conversions in the PLL by 
accounting for the transfer function H(P), where P = ^ is the Laplace operator. If 
H(P) is a low pass filtering function with cut-off frequency f c , the frequency at the 
output of the mixer + filter stage is such that 

A* = /b(*)//o = QiW ~ Pi/Q*\ ^ /c//o, Pi and qi integers. (6) 
The beat frequency results from the continued fraction expansion of the input 

frequency ratio 

v = f(t)/fo = [a ;ai,a 2 , . . .a i; a, . . .], (7) 

where the brackets mean expansions ao + 1 / (ai + 1 / (a 2 + l / . . .+l/(oj+l/(a+. • •))))• The 
truncation at the integer part a = [^r] defines the edges of the basin for the resonance 
Pi/qf, they are located at v x = [a ; ai, a 2 , . . . , a iy a] and u 2 = [a ; cti, a 2 , . . . a.j — 1, 1, a] pQ. 
The two expansions in v\ and u 2 , prior to the last filtering partial quotient a, are the 
two allowed ones for a rational number. The convergents Pi/qi at level i are obtained 
using the matrix products 

(8) 

= ^ / 7~v\ ^~ i so that 

qi{2a+l)-q a ' 

the width of the basin 



a 


1 " 




d\ 


1 " 




^ 1 




Pi 


Pi-l 


1 


_ 




1 


_ 




10 






Qi-i _ 



Using (jHJ), one can get the fractions V\ and u 2 as v x 
with the relation relating convergents (ptqi-i —pi-iqi 



^ and v 2 

la Z 

= f-iV- 1 , 



of index i is \v\ — v 2 \ 



2a+l 



whenever a > 1. 



<3a(<?a+(2a+l)qi) q a qi 

In previous publications of one of the authors a phenomenological model for 
1/f noise in the PLL was proposed, based on an arithmetical function which is a 
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logarithmic coding for prime numbers P0,[2]- If one accepts a coupling coefficient 
evolving discontinuously versus the time n as K = K A(n), with A(n) the Mangoldt 
function which is ln(p) if n is the power of a prime number p and otherwise, then the 
average coupling coefficient is K Q and there is an arithmetical fluctuation e(t) 



V>(f) = £A(n)=f(l + e(f)), 



n=l 



te(t) = -ln(27r) - - ln(l - t~ 2 ) 



t p 

p p 



(9) 



The three terms at the right hand side of te(t) come from the singularities of the Riemann 
zeta function ((s), that are the pole at s — 1, the trivial zeros at s = —21, I integer, and 
the zeros on the critical line 3?(s) = |[T|. Moreover the power spectral density roughly 
shows a 1/ f dependance versus the Fourier frequency /. This is the proposed relation 
between Riemann zeros (the still unproved Riemann hypothesis is that all zeros should 
lie on the critical line) and 1/f noise. 

We improved the model by replacing the Mangoldt function by its modified form 
b(n) = A(n)(p(n)/n, with <p{n) the Euler (totient) function [HI]. This seemingly 
insignificant change was introduced by Hardy |Hj in the context of Ramanujan sums 
for the Goldbach conjecture and resurrected by Gadiyar and Padma in their recent 
analysis of the distribution of pairs of prime numbers [HZ] , Then by defining the error 
term e#(t) from the cumulative modified Mangoldt function 



B(t) = £&(T0=f(l + e B (t)), 



(10) 



n=l 



its power spectral density SbU) — exhibits a slope close to the Golden ratio 
(VE - l)/2 ~ 0.618 (see Fig. 1). 
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Figure 1. Power spectral density of the error term in the modified Mangoldt function 
b(n) in comparison to the power law 1/ f 2a , with the Golden ratio a = (y/5 — l)/2. 
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The modified Mangoldt function occurs in a natural way from the logarithmic 
derivative of the following quotient 

7 i~\ _ <M _ y MO nil 

since — = X)n>i ^r- This replaces the similar relation from the Riemann zeta 

function where = E„>i ^r- 

In the studies of 1// noise, the fast Fourier transform (FFT)plays a central role. 
But the FFT refers to the fast calculation of the discrete Fourier transform (DFT) with 
a finite period q = 2 l , I a positive integer. In the DFT one starts with all q th roots of 
the unity exp(2i7cp/q), p = 1 . . . q and the signal analysis of the arithmetical sequence 
x(n) is performed by projecting onto the nth powers (or characters of Z/qZ) with well 
known formulas. 

The signal analysis based on the DFT is not well suited to aperiodic sequences with 
many resonances (naturally a resonance is a primitive root of the unity: (p, q) = 1), and 
the FFT may fail to discover the underlying structure in the spectrum. We recently 
introduced a new method based on the Ramanujan sums defined in (J3J) |lUj . 

Mangoldt function is related to Mobius function thanks to the Ramanujan sums 
expansion found by Hardy 



Kn) = M0 A(n) = f ; ^M Ci(n) . (12) 



n 



We call such a type of Fourier expansion a Ramanujan- Fourier transform (RFT). General 
formulas are given in our recent publication and in the paper by Gadiyar ^2] . This 
author also reports on a stimulating conjecture relating the autocorrelation function of 
b(n) and the problem of pairs of prime numbers. In the special case fll2j) . it is clear that 
/i(q , )/0(g) is the RFT of the modified Mangoldt sequence b{n). 

Using Ramanujan- Fourier analysis the l/f 2a power spectrum gets replaced by a 
new signature shown on Fig. 2, not very different of //(g) / '<fi(q) (up to a scaling factor). 

2.2. The hyperbolic geometry of phase noise and 1/f frequency noise 

The whole theory can be justified by studying the noise in the half plane H = {z = 
v + iy, i 2 = —1%, y > 0} of coordinates v = -j- and y — ^ > and by introducing 
the modular transformations 

/ \ / Vi Z "I" Pi „ / / 1 /-in\ 

z y(z) = z = — , p i q i - pm = 1. (13) 

QiZ + Qi 

The set of images of the filtering line z = v + % under all modular transformations can 
be written as 

^~{t + tL)\ = T2- ( 14 ) 



% The imaginary symbol i should not be confused with the index i in integers pi, qi and in related 
integers. 
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Figure 2. Ramanujan-Fourier transform (RFT) of the error term (upper curve) of 
modified Mangoldt function b(n) in comparison to the function fi(q)/cf>(q) (lower curve). 



Equation (|T3|) defines Ford circles (see Fig. 3) centered at points z = ^ + ^ with 
radius ^ CH]- To each — a Ford circle in the upper half plane can be attached, which 
is tangent to the real axis at u — Ei . Ford circles never intersect: they are tangent 
to each other if and only if they belong to fractions which are adjacent in the Farey 




q 1/4 1/3 1/2 2/3 3/4 . 

Figure 3. Ford circles: the mapping of the filtering line under modular 
transformations (|13|l . The arrows indicates that Ford circles were used as an integration 
path by Rademacher to compute the partition function p(n) . 

hyperbolic geometry. A basic fact about the modular transformations (fl3|) is that they 
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form a discontinuous group T ~ SL(2, Z) /{±i }, which is called the modular group. The 
action of T on the half-plane H looks like the one generated by two independent linear 
translations on the Euclidean plane, which is equivalent to a tesselation the complex 
plane C with congruent parallelograms. One introduces the fundamental domain of V 
(or modular surface) F = {z E H : |z| > i, |v| < 4}, and the family of domains 
{7(F), 7 E r} induces a tesselation of H [14J. 

It can be shown ^],[Hj that the noise amplitude is a particular type of solution 
of the eigenvalue problem with the non- Euclidean Laplacian A = y 2 (-^2 + J^)- The 
solution corresponds to the scattering of waves in the fundamental domain F. It can be 
approximated as a superposition of three contributions. The first one is an horizontal 
wave and is of the power law form y s , the second one is also horizontal wave of the form 
S(s)y 1 ~ s and corresponds to a reflected wave with a scattering coefficient of modulus 
\S(s)\ = 1, whereas the remaining part T(y,u) is a complex superposition of waves 
depending of y and the harmonics of exp(2i7rz/), but going to zero for y — ► 00. Extracting 
the smooth part in S(s) one is left with a random factor which is precisely equal to the 
function Z(2s — 1) defined in (fTT|) as the quotient of two Riemann zeta functions at 
2s — 1 and 2s, respectively. An interesting case is when s is on the critical line, i.e. 
s = \ + ik in which case the superposition T(y, v) vanishes and the reflexion coefficient 
is 

S(k) = exp[2i6(k)], with 6'{k) = at s = i + ik. (15) 

The scattering of waves from the modular surface is thus similar to the phase-locking 
model plotted in the set of equations (J5|)- (fTTj) . It explains the relationship between the 
hyperbolic phase and the 1/ / noise found in the counting function 6'{k). The phase 
factor 6{k) is represented in Fig. 4. 

3. Quantum phase-locking 

3.1. The quantum phase operators 

Going back to the quantum definition of phase states announced in the introduction 
one calculates the projection operator over the subset of phase-locked quantum states 
K) as 

^ ock = E KW P \ = -E c > - OlnX'l, (16) 

V 1 n,l 

where the range of values of n,l is from to (p(q). Thus the matrix elements of the 
projection are q(n\P q \l) = c q (n — I). This sheds light on the equivalence between 
cyclotomic lattices of algebraic number theory and the quantum theory of phase-locked 
states. 
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Figure 4. The phase angle 9(k) for the scattering of noise waves on the modular 
surface. Plain lines: Exact phase factor. Dotted lines: Approximation based on the 
quotient of two Riemann zeta functions^ . 



The projection operator over the subset of pairs of phase-locked quantum states 
\9' ) is calculated as 



1 



„ - ' ' ' (17) 

P,P ^ n,l 

where the notation p,p means that the summation is applied to such pairs of states 
satisfying (£Q). The matrix elements of the projection are q(n\P^ aas \l) = k q (n,l), which 
are in the form of so-called Kloosterman sums jTH] 



1m 

k g (n,l) = E ex P[ {pn-pl)}. 

ViV ^ 



Kloosterman sums k q (n,l) as well as Ramanujan sums c q (n — I) are relative integers. 
They are given below for the Hilbert dimensions q = 5(0(5) = 4) and q = 6(0(6) = 2). 
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One defines the quantum phase-locking operator as 



2^. 



(19) 
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The Pegg and Barnett operator [7] is obtained by removing the coprimality condition. 
It is Hermitian with eigenvalues 6 p . Using the number operator N q = X)n=o n l n )( n l a 
generalization of Dirac's commutator [Q q , N q ] = —i has been obtained. 

Similarly one defines the quantum phase operator for Kloosterman pairs as 

Qpairs = £ ^ ^| = ^ppairs ^ ^ = ^ P (2Q) 
p,p 1 

The phase number commutator for phase-locked states calculated from ()19j) is 

c io Ck = [e io Ck) Nq] = 1J2(1- n)c q (n - l)\n)(l\, (21) 

1 n,l 

with antisymmetric matrix elements (/|C^ ock |n) = -(/ — n)c q {n — I). 

For pairs of phase-locked states an antisymmetric commutator CP airs similar to (j21j) 
is obtained with k q (n, I) in place of c q {n — I). 

3.2. Phase expectation value and variance 

The finite quantum mechanical rules are encoded in the expectation values of the phase 
operator and phase variance. 




Figure 5. Oscillations in the expectation value (|23|l of the locked phase at j3 = 1 
(dotted line) and their squeezing at /? — (plain line). The brokenhearted line which 
touches the horizontal axis is wA(q)/hiq. 



Rephrasing Pegg and Barnett, let us consider a pure phase state |/) = J2n=o u n\ n ) 
having u n of the form 

u n = (l/v/g)exp(m/3), (22) 
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where (3 is a real phase parameter. One defines the phase probability distribution 
(0' p \f) 2 , the phase expectation value (6g° ck ) = J2 p Qp(Q p \f) 2 , and the phase variance 

(AQ^lock = Ep(dp _ (e lock ))2( ^| /} 2 Qne gets 

( @ i° Ck ) = - 2 Y. C S - n)exp[iP(n ~ I)), (23) 

" n,l 

(AQ 2 Q ) lock = 4(6' ock ) + ^^«e 9 ) - 2tt), (24) 

with the modified expectation value (0g ) = ^J2n,iC q (l — n) exp[if3(n — I)}, and the 
modified Ramanujan sums c q {n) = Y^pip/l) 2 exp(2z7rm 2 ). 

Fig. 5 illustrates the phase expectation value versus the dimension q for two 
different values of the phase parameter (3. For (3=1 they are peaks at dimensions 
q = p r which are powers of a prime number p. The most significant peaks are fitted by 
the function 7rA(g)/ lng, where A(g) is the Mangoldt function introduced in Q of Sect. 2. 
This observation provides the link between the arithmetical hyperbolic viewpoint and 
the quantum one. A deepest explanation based on the relation with quantum statistical 
mechanics and the work of Bost and Connes can be found in ^7j. For (3 = the peaks 
are smoothed out due to the averaging over the Ramanujan sums matrix. Fig. 6 shows 
the phase expectation value versus the phase parameter (3. For the case of the prime 
number q = 13, the mean value is high with absorption like lines at isolated values of (3. 
For the case of the dimension q = 15 which is not a prime power the phase expectation 
is much lower in value and much more random. 

Fig. 7 illustrates the phase variance versus the dimension q. Again the case 
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(3 = 1 leads to peaks at prime powers. Like the expectation value in Fig. 5, it is 
thus reminiscent of the Mangoldt function. Mango ldt function A(n) is defined as Inp 
if n is the power of a prime number p and otherwise. It arises in the frame of prime 
number theory pQ from the logarithmic derivative of the Riemann zeta function ((s) 
as — = Y^=o ^r- average value oscillates about 1 with an error term which is 
explicitely related to the positions of zeros of £(s) on the critical line s = |. The error 
term shows a power spectral density close to that of 1/ f noise p. It is stimulating to 
recover results reminding prime number theory in the new context of quantum phase- 
locking. 

Finally, the phase variance is considerably smoothed out for (3 = it and is much 
lower than the classical limit tt 2 /3. The parameter (3 can thus be interpreted as a 
squeezing parameter since it allows to define quantum phase-locked states having weak 
phase variance for a whole range of dimensions. 




_0,5 i i i i i i i i i i I 

5 10 15 20 25 30 35 40 45 50 
dimension 

Figure 7. Phase variance versus the dimension q of the Hilbert space. Plain lines: 
(3=1. Dotted lines: /3 = w. 



3.3. Towards discrete phase entanglement 

The expectation value of quantum phase states can be rewritten using the projection 
operator of individual phase states 7i p = \0' p ){6' \ as follows 

<er> = EUW W> = EWKlf)- (25) 
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dimension q 



Figure 8. Phase expectation value versus the dimension q for pairs of phase- locked 
states. Plain lines: (3 — 0. Dotted lines: (3 = 1. 

This suggests a definition of expectation values for pairs based on the product 
follows 

(6f ) = 5>p</lvr*l/>- (26) 

p,p 

It is inspired by the quantum calculation of correlations in Bell's theorem [THI- Using 
pure phase states as in (J22j) we get 

(©f rs ) = % E K(n, 1) ex P [z/3(n - I)}, (27) 
where we introduced generalized Kloosterman sums 

k q (n,l) = Y^p,pP ex P[—(p — P)(l — n )}- These sums are in general complex numbers 
(and are not Gaussian integers). The expectation value is real as expected. In Fig. 
8 it is represented versus the dimension q for two different values, (3 = and (3=1, 
respectively. Note that the pair correlation ()26|) is very strongly dependent on q and 
becomes quite huge at some values. 

This result suggests that a detailed study of Bell's type inequalities based on 
quantum phase-locked states, and their relationship to the properties of numbers, should 
be undertaken. Calculations involving fully entangled states 

|/) = -£|0 P) l>®|%,2>, (28) 
" p,p 

have to be carried out. This is left for future work. 
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Table 1. (Z/7Z)* is a cyclic group of order (f)(7) = 6. 
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Table 2. (Z/3 2 Z)*, is a cyclic group of order (f)(9) = 6. 
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Table 3. (Z/8Z)* lias a largest cyclic group of order A(8) = 2. 
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1 



S.Jf.. The discrete phase: cycles in Z/qZ 

There is a scalar viewpoint for the above approach, which emphasizes well the intricate 
order of the group Z/qZ, the group of integers modulo q. One asks the question: what 
is the largest cycle in that group. For that purpose one looks at the primitive roots, 
which are the solutions g of the equation 

g a = l(mod q), (29) 

such that the equation is wrong for any 1 < a < q — 1 and true only for a = q — 1. 
If q=p, a prime number, and p = 7, the largest period is thus 0(jo)=p-l=6, and the 
cycle is as given in Table I. If q = 2, 4, q = p r , a power a prime number > 2, or 
q = 2p r , twice the power of a prime number > 2, then a primitive root exists, and the 
largest cycle in the group is 4>(q). For example g = 2 and q = 3 2 leads to the period 
0(9) = 6<g— 1 = 8, as shown in Table II. Otherwise there is no primitive root. The 
period of the largest cycle in Z/qZ can still be calculated and is called the Carmichael 
Lambda function X(q). It is shown in Table III for the case g = 3 and q — 8. It is 
A(8) = 2 < (f)(8) = 4<8 — 1 = 7. Fig. 9 shows the properly normalized period for the 
cycles in Z/qZ. Its fractal character can be appreciated by looking at the corresponding 
power spectral density shown in Fig. 10. It has the form of a 1/ f a noise, with a = 0.70. 
For a more refined link between primitive roots g, cyclotomy and Ramanujan sums see 
also [HI. 

4. Conclusion 

In conclusion, we explained how useful could be the concepts of prime number theory 
in explaining various features of phase-locking at the classical and quantum level. In 
the classical realm we reminded the hyperbolic geometry of phase, which occurs when 
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Figure 10. FFT of the normalized Carmichael lambda function. The staight line has 
slope -0.70. 

one accounts for all harmonics in the mixing and low-pass filtering process, how 1/f 
frequency noise is produced and how it is related to Mangoldt function, and thus to the 
critical zeros of Riemann zeta function. Then we studied several properties resulting 
from introducing phase-locking in Pegg-Barnett quantum phase formalism. The idea of 
quantum teleportation was initially formulated by Bennett et al in finite-dimensional 
Hilbert space jTHj, but, yet independently of this, one can conjecture that cyclotomic 
aspects in phase-locking could play an important role in many fundamental tests of 
quantum mechanics related to quantum entanglement. Munro and Milburn [20] already 
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conjectured that the best way to see the quantum nature of correlations in entangled 
states is through the measurement of the observable canonically conjugate to photon 
number, i.e. the quantum phase. In their paper dealing with the Greenberger-Horne- 
Zeilinger quantum correlations, they presented a homodyne scheme requiring discrete 
phase measurement. We expect that the interplay between quantum mechanics and 
number theory will appear repetitively in the coming attempts to manipulate quantum 
information [2~T] . 
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